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considering the different operational modes of the instrument. 



Keywords Polarization, Optical • Instrumental Effects • Magnetic Fields, 
Photosphere 



J.M. Borrero, S. Tomczyk, A. Norton*, T. Darnell 

High Altitude Observatory, National Corporation for 

Atmospheric Research, 3450 Mitchell Lane, Boulder CO 80301, USA 

E-mail: borrero@ucar.edu, tomczyk@ucar.edu, norton@ucar.edu, tdar- 

nell@ucar.edu 



J. Schou, P. Scherrer, R. Bush, Y. Liu 

Hansen Experimental Physics Laboratory, Stanford University 
Stanford, CA 94305, USA 

E-mail: schou@stanford.edu, pscherrer@stanford.edu, rock@stanford.edu, 
yliu@stanford.edu 



* Current address: National Solar Observatory, 950 N Cherry Avenue, Tucson AZ 
85719 



2 



1 Introduction 

The Solar Dynamics Observatory is scheduled to launch from Cape Canaveral 
in September 2008 on an Atlas V Booster. Onboard this satellite there will 
be several instruments dedicated to the study of the solar photosphcric and 
coronal magnetic fields and their relation to the interplanetary medium, space 
weather and terrestrial climatology. The Helioseismic and Magnetic Imager 
is an instrument developed by a collaboration between Stanford University, 
Lockheed-Martin Solar and Astrophysics Laboratory and the High Altitude 
Observatory. HMI consists of a combination of a Lyot filter and two Michelson 
interferometers that will measure the full Stokes vector I = {I,Q,U,V) at 
several wavelengths across the Fe I 6173.3 A line. HMI will have twin CCD 
cameras that will operate independently. The first one, to be referred to as 
the "Doppler camera" , will be devoted to the measurement of the line-of- 
sight component of the magnetic field and velocity, like MDI on SOHO. The 
second camera, hereafter referred to as the "Vector Camera", will measure 
the vector magnetic field and line-of-sight velocities. Each CCD consists of a 
4096 X 4096 pixel array with a pixel size of 0.5 arcsec. Together with SDO's 
geosynchronous orbit, this will allow for a near continuous monitoring of the 
solar (full disk) vector magnetic field. 

In this paper we study the accuracy to which we can retrieve the mag- 
netic field vector and other physical parameters of the solar photosphere 
using the data from the HMI vector camera. Details about how the vector 
camera records the data are presented in Section 2. In Section 3 we describe 
the main sources of uncertainty in the data (photon noise, limited spectral 
resolution, solar oscillatory modes) and how we simulate HMI data by taking 
into account those effects. Once we account for those processes, we are able 
to use a Stokes synthesis code to compute polarization profiles of the Fe I 
6173.3 A line as if they were observed by HMI. A large database of profiles is 
built by using an array of atmospheric parameters. Using a Stokes inversion 
code (described in Section 4) the physical parameters used in the synthesis 
can be recovered. The errors in the retrieval depending upon the choice of the 
polarization modulation are studied in Section 5. The effects of the photon 
noise are considered in Section 6. Section 7 investigates the improvement in 
the retrieval when the data are averaged in time. Our conclusions are finally 
given in Section 8. 

2 HMI Operational Modes 

Line selection: The Helioseismic and Magnetic Imager will obtain two-dimensional 
images of the full solar disk at several wavelengths across the neutral iron 
line at 6173.3 A. This spectral line was chosen after a careful comparison 
with other candidates: Ni I 6767.8 and Fe I 6302.5 (see Norton et al. 2006). 
This line possesses a Lande factor of ge// =2.5 and therefore is more suit- 
able for the investigation of the magnetic field than its MDI counterpart Ni I 
(ge// =1.5) 6767.8 A (see Scherrer et al. 1995). In addition, its large slope in 
the wings ensures a high accuracy in the derived velocity signals (see Cabrera 
Solana, Bellot Rubio, and del Toro Iniesta, 2005). It is also relatively free of 
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blends which makes this hne very convenient for HMI purposes, since the 
geosynchronous orbit of SDO will induce Doppler shifts as large as ±4 km 
s~^ due to spacecraft velocities. 

Filter profiles: The combined Lyot-Michelson filter system in HMI produces 
a transmission profile with a FWHM of 76 mA. Tuning positions are sepa- 
rated by 69 mA. Although there are a number of different configurations and 
possible tuning positions available, in this work, we will restrict ourselves 
to studying a configuration in which the spectral line is scanned across five 
wavelength positions with the middle one located at the central laboratory 
wavelength Aq of the spectral line. Figure 1 displays the filter profiles at each 
wavelength position. In the remainder of the paper we will denote each of 
the filter profiles as: t/i(A), where i — 1, M — 5. 

Filter sequence and polarization modulation: For each tuning position (Ai) 
and time {tn), HMI will measure a linear combination [0(Ai, of the com- 
ponents of the solar Stokes vectoiQ: 

4 

Op(A„i„) = ;^Mpfe7f '^(A„i„) (1) 

k=l 

Where Mpk are the elements of the M modulation matrix with p = 
1,...,P. A4 is therefore a P x 4 matrix. Before shifting to the next tuning 
position (i + 1), HMI will measure the next linear combination of Stokes 
profiles, as given by the modulation matrix, at the same wavelength: 

4 

Op+i{Xi,tn+l) — ''^^Mp-\-l^kIk^^{K,tn+l) (2) 
k=l 

Once this is done the wavelength position is shifted to i -I- 1 and the pro- 
cess for p,p+l is repeated until i — M. The observation sequence continues 
by changing Equations (1) - (2) to p + 2,p + 3 and so forth until the last scan 
is made with P —1, P. Throughout this process, the time index has run from 
n = 1, A/ X P. 

I°^^{Xi,tn) results from applying the HMI filter function at the tuning po- 
sition Ai, Gi{Xj) to the original {i.e.: solar) Stokes vector I°^^{Xj,tn) sampled 
in a continuous wavelength grid Xj at a time t„. 

r'^^(A„ i„) = I r^^(A„ i„)e^(Aj - AOdA, (3) 

In order to obtain the observed Stokes vector at each wavelength position, 
one must first demodulate the vector 0{Xi) from Equation (1). 



^ In compact notation, the components of the Stokes vector are expressed as Ik 
{k = 1, ...4) with: h = l,h = Q, h = U,h = V 
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Fe I 6173.35 A 
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Fig. 1 HMI filter functions Gi{Xj) at five different tuning positions. An example 
of the intensity profile for the line Fe I 6173.3 A is given in solid. 

^'"'{X^) = M-'0{X,) (4) 
or equivalently, as follows. 

p 

ifW) ^Y.^kpO^^>^M (5) 

Note that both the vector 0{\i) and the observed Stokes vector I°'^*'(Ai) 
(or any of its k components) are gathered over APM seconds since each 
element Op was acquired at a diiferent time i„. Therefore, what is observed 
by HMI is some sort of temporal average of the solar Stokes vector. 

A number of possible filter sequences and polarization modulations have 
been suggested. In this work we will restrict ourselves to consider only three of 
thcrr0. They will be referred to as: modulations A, B or C. Sample framelists 
for them are given in Tables I, II, and III. Referring to those Tables and Equa- 
tion (1) it is possible to construct the modulation matrices M. corresponding 
to modulation schemes A, B, and C. 

^ All these modulation schemes take into account that the Doppler camera will 
scan /± V in wavelength with a cadence of 50 seconds or better. For the vector cam- 
era it would have been more desirable to measure all of polarization states at each 
tuning position (Ai) before shifting to the next wavelength. However, this would 
require too many motor rotations between observations which may risk mechanical 
failure. 
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Table I Pramelist sample for modulation scheme A. The linear combination of 

the components of the Stokes vector measured at each time are given by: Oi = 
I + aQ + bV, O2 = I - aQ + bV, O3 = I + aU - bV, O4 = I - aU - bV. Constant 
values for the linear combinations are: a = 0.81 and b = 0.588. Time is in units 
of At = four seconds. At is determined by CCD readout time, the time needed 

to select a polarization state and to shift the filter position. Note that the total 
time needed for modulation A or B to measure the full Stokes vector across five 
wavelength points is PM At = 20At = 80 scchikIs. 
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(6) 



(7) 



(8) 



It is important to mention that modulations A and B require 80 seconds 
to record the full Stokes vector at all tuning positions, while modulation 
C needs 120 seconds. A and B demodulate according to Equations (4) - (5) 
since A4 is a 4x4 square matrix (sec Equations (6) -(7)) and its inverse is 
calculated easily. On the contrary, modulation scheme C requires some other 
consideration since Mc (Equation (8)) is not square. This is due to the 
fact that modulation C gathers more information than actually needed. In 
particular. Stokes /(A) is obtained three times. For simplicity, the average 
/(A) will be used. 

It is also important to note that all schemes have a polarimetric efficiency 
of e = 1/a/3 when measuring Stokes Q, U, and V. According to del Toro 
Iniesta and CoUados (2000) this translates into an error (due to photon noise) 
ar. 



1 

2 



(9) 
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Table II Same as Table I but for modulation scheme B. 
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Table III Framelist sample for modulation scheme C. The linear combination of 
the components of the Stokes vector measured each time are given by: Oi = I + Q, 
O2 = I - Q, Os = I + U, O4 = I -U, Oi = I +V, Oe = I - V. Note that the 
total time needed by modulation C to measure the full Stokes vector across five 
wavelength points is PMAt — 30At = 120 seconds. 
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where cr^ is roughly the noise level in Op (defined in Equation (1)). Note that 
scheme C, vi^here P = 6, introduces a factor 2/3 less noise than schemes A 
and B (P = 4). However, scheme C takes a factor 3/2 longer to be completed 
(120 seconds instead of 80 seconds) so that the noise in the Stokes parameters 
ai per unit time is the same for all three schemes. 



3 Simulated HMI Observed Profiles 

In order to simulate Stokes profiles as if they were observed by the HMI vector 
camera, a set of parameters [Xr(t„)] describing the conditions present in the 
solar photosphere at a time tn, is chosen. In the following we will assume that 
the solar photosphere is described by a Milne-Eddington (ME) atmosphere 
(i.e.: physical parameters arc constant with depth in the photosphere: see 
Landi Degl'Innocenti (1991) for details). This set of parameters is presented 
in Table IV. Xr(t„) is used to solve the radiative transfer equation (RTE) 
according to the ME analytical solution given by Landi Degl'Innocenti (1991) 
or Del Tore Iniesta (2003). This provides the Stokes vector emerging from a 
solar magnetized atmosphere [I™8(Xr)]. In addition, a contribution from a 
non-magnetic atmosphere (I"™) is added using the filling factor a (fractional 
area of the resolution element occupied by the magnetic atmosphere) in order 
to produce the simulated observed profiles: 

r^^(A,-,t„) = aI'"«(A„Xr[f„]) + (1 - a)r™(A,) (10) 
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Table IV Set of Milne-Eddington physical parameters Xr used to produce syn- 
thetic profiles. 



Physical Parameter 
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or cquivalcntly, taking into account that the non-magnetic atmosphere does 
not produce polarization signals: 



/jobs \ 




/jmg 




— a 




^obs 


[/mg 









+ (1 




(11) 



1° (Aj, t„) is sampled using a grid of 300 points in wavelength (as an 
approximation to a continuous function; see Equation (3)). This Stokes vector 
enters into Equation (3), where the HMI filter profiles Qi are applied and 
l°^^{\i,tri) (sampled with a coarse grid of five wavelength points) is obtained. 
According to the selected polarization modulation, M, a linear combination 
of the components of the Stokes vector, Op{\i, fn), is produced (see Equation 
(I))- 



3.1 PHOTON NOISE 

Considering that the full-well depth of the CCD camera is approximately 
100 000 detected photons, it is reasonable to limit the exposure time in order 
to reach 4/5 of that value in order to avoid saturation effects. Photon noise 
a is added into Op(Ai,t„) considering that 4/5 of the CCD full well depth 
is achieved for the brightest solar feature (i.e.: quiet Sun at disk center). 
For profiles intended to simulate other solar features, e.g. sunspots. photon 
noise is added according to the expected intensity level. The noise is added 
using random number generation following Poisson's statistics. This is often 
referred to as sh,ot noise. 

Note that this noise level a translates into different noise for the po- 
larization signals ai (see Equation (9)). Since we add noise into Op{\i), Ui 
propagates naturally after demodulating (Equations (4) -(5)). 



3.2 SOLAR P MODES 



The readout time of the 4096 x 4096 CCD camera is about 2.7 seconds. To- 
gether with the time required to rotate the waveplates and some idle time, the 
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total time lapse before Op+i{Xi,tn + At) is measured is: At = four seconds. In 
this time interval the physical conditions in the solar photosphere might have 
changed: Xr(i„ + At) ^ Xr(i„), and therefore the observed l°^^{Xj,tn + At) 
is likely to be different from I°^^(Aj, 

There are many processes responsible for the variation of the physical 
conditions of the solar photosphere on time scales of seconds. In this paper 
we will focus of the effects of the p modes only. This means that from the set 
of parameters Xr(f„) = [B,7,V',v^'^^,S,S], only v^'^® will be considered to 
change in time. In particular: 

v^°'(t)=vLOS_ + ^(t) (12) 

where Vpi^^^jj^g^ corresponds to a real plasma velocity, while £,{t) is the velocity 
associated with the solar p modes. An example of the function ^{t) is pre- 
sented in Figure 2 (top panel). This was obtained from a two hour average 
of quiet-Sun MDI data. The rms amplitude of the oscillation is about 250 m 

Long time-scale p-mode effect A close inspection of Figure 2 (middle panel) 
indicates that, for example using modulation scheme C, during the time 
interval [to, to + 2MAt] the HMI vector camera will measure: Co (A) = /(A) + 
(5(A) and Oi(A) = /(A) — Q(A). This allows an acquisition of both /(A) and 
(5(A). During the next time interval [to + 2MAt,to + AM At] HMI records: 
e>2(A) = I{X) + U{X) and 03(A) = I{X)-U{X). This provides /(A) and U{X). 
Due to the temporal evolution of the solar j>mode velocity [£,{t)] during 
this time, Q{X) and U{X) are affected, on average, by different velocities 
and therefore they will be shifted in wavelength with respect to each other. 
Hereafter this effect will be referred to as the long time scale p-mode effect. 

Short time scale p-mode effect Another effect associated with the solar p 
modes can be visualized by looking at Figure 2 (bottom panel). The time 
to is when the vector camera records Op{Xi). Immediately after: to + i, is 
when HMI measures the state Op+i{Xi): Let us consider for example that 
modulation scheme C is used. In this case t = At = four seconds. We can 
thus write: 

02p-iiM,to) = lf'{X^,to) + i;X\{k,to) m= 1,2,3 (13) 
and a time t later the orthogonal state is measured: 

02p {Xi ,to + t) = If' {Xi ,to + t)- i;X\ {Xi ,to + t) (14) 

In order to isolate the polarization signal one would normally proceed 
as follows: 

02p-i{Xi, to) — 02p{Xi,to +i) — I°^^{Xi,to) — I°^^{Xi,to + i) ^^^^ 

+ /p+i(Aj,fo) + Ip+\{\,to -\-i) 
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Fig. 2 Top panel: dashed line describes the function ^(t) (Equation (12)) represen- 
tative of the solar p modes as measured by the MDI high resolution field of view. 
The vertical solid line indicates a randomly chosen time to at which HMI starts 
acquiring data. The solid line indicates the variation of the solar p-mode velocity 
during the time interval [to, to + PMAt] . This is the time required to record the 
full Stokes vector across M tuning positions. Middle panel: close up of the function 
^(t) over the data acquisition time interval (solid line in the upper panel). The 
vertical dashed lines separate the time intervals in which Op (A) and Op+i[\) are 
recorded at all tuning positions, taking a total of [to, to + 2MAt] indicated by the 
solid line. Bottom panel: close up of the solid line in the middle panel. Each vertical 
dashed line separates the time interval in which Op and Op+i are taken at each 
tuning position (Aj). Between Op and Op+i, there is a time delay At. 



Both Ii and /,„ have changed between to and to + i. The intensity profiles 
{i.e.: Stokes / or Ii) have a much larger amphtude that the polarization 
profiles /„. Therefore we can ignore changes in and consider that the 
solar p modes have shifted only the intensity profiles: 



02p-i{Xi, to) - 02p{Xi, to + i) ~ 2I^\\{Xi,to) + 



dt 



i (16) 



t=to 



xtalk I\ 



where the term indicated as xtalk in Equation(16) will be referred to as short 
time-scale p-mode effect. An order-of-magnitude estimation of this effect can 
be made by modeling the intensity profile as a Gaussian and the temporal 
evolution of the solar p modes in the following form: 
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lf%X,t) = 1-Ic exp (-/J[A - Ao(l - ^)]^) (17) 

^{t) = Vocos{ — ) (18) 

We can easily calculate the cross talk from Stokes / into the polarization 
profiles as: 



A-Ao(l-^) 
c 



Using the following typical values for the neutral iron line Fe I 6173.3 

A taken from Graham et al. (2002, 2003) and Norton et al. (2006): Ic = 0.5, 
Vo = 250 m s-i (see Figure 2; top panel), (3 = 385 A'^, and T = 300 s, 
we have evaluated Equation (19) and plotted it in Figure 3 in units of per- 
centage of the continuum intensity. As expected, the maximum cross talk is 
produced in the line wings because that is where changes in Ji are maximum 
when the velocity changes. The maximum cross talk is produced when the 
p-mode velocity changes sign. Figure 3 also shows that along the time axis 
the cross talk changes sign for a particular wavelength {e.g.: in the increas- 
ing region of the p-mode velocity the cross talk has opposite sign than in the 
decreasing velocity region), implying that integrating over time can reduce 
the cross talk level (sec Section 6J. The gray scale in Figure 3 is set so that, 
in the best case scenario, where t = At = four seconds (modulation scheme 
C), we have an estimated cross talk of ~ 0.15 x 4 = 0.6 %, while in other 
cases, where t ~ 50 — 80 seconds, the amount of cross talk could be as much 
as 7-9 %. The modulation matrix M. in Equations (6) -(8) defines how far 
away in time the measurements Op (that will be combined to produce each 
Stokes profile) are carried out and therefore defines i. 



Besides this analytical estimate, we have carried out a more detailed 
analysis by taking into account all possible effects and a realistic p-mode 
function £^{t) (Figure 2 top panel). For this analysis we used 20 000 Stokes 
profiles P'"^(Aj,t„) that were converted into P''^(Aj) according to Section 
2. Polarization levels, defined as max|y| and max| + C/^| in units of 
percentage of the continuum intensity, are compared before and after. Results 
are displayed in Figure 4. The horizontal axis refers to the polarization levels 
from the orig inal profiles: I°'^"(A,,t„) as in Equation (3). The vertical axis 
represents the polarization levels of the profiles once the effects of the solar 
p modes have been taken into account: I°^''(Ai) as in Equation (4). Random 
initial times arc used for sampling the function ^(f) (sec Figmc 2: top 
panel). The experiment is repeated for all modulation schemes (see Section 
2). As can be seen, modulation scheme A introduces a large amount of cross 
talk into the circular polarization while keeping the cross talk into the linear 
polarization within reasonable limits. 
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Fig. 3 Contour plot of the amount of cross talk per unit time (in units of % of the 
continuum intensity) induced into the polarization profiles from the total intensity 
as a consequence of the solar p modes (see Equation (19) and text for details). 
Horizontal axis indicates relative position with respect to the central laboratory 
wavelength. Vertical axis indicates time in which the measurement is taken within 
one oscillatory period of 300 seconds. Relative amplitudes are indicated by the gray 
scale bar at the top of the plot. For visualization purposes, a magnified example 
of the Fe I 6173 A spectral line has been over plotted in dashed lines. Temporal 
evolution of the p-mode velocity ^{t) (see Equation (18)) is indicated by the solid 
lines. 



As already mentioned, this effect is ultimately related to the modulation 
matrix M- In scheme A, Ci(Aj, to) = I — 0.81Q + O.SSSV it is the first mea- 
surement. The next one corresponds to 02{Xi,to + At) = / — 0.81(3 + 0.588^ 
and therefore. Stokes Q{Xi) can be constructed with measurements done only 
At = four seconds apart from each other, thereby minimizing the amount of 
cross talk from Stokes /. However, Stokes V is built out of a combination of 
measurements taken far away in time. For example: C'i(Aj, to) = + 0.81Q + 
0.588y is measured at to, but Oz{\i, to + [2M + l]At) = 1 + 0.81C/ - 0.588^ 
is recorded 44 seconds later (see Table I). Dming this time the solar p modes 
might introduce large levels of cross talk from Stokes / into Stokes V. This 
also explains why scheme B introduces large cross talk into the linear po- 
larization while leaving the circular polarization almost unchanged. As ex- 
pected, Modulation scheme C performs well both in linear and circular po- 
larization. The obvious drawback to Scheme C is that it takes longer to scan 
across the line: 6MAt seconds as compared to AM At seconds from schemes 
A and B. Thus, scheme C minimizes the short time scale p-mode effect but 
maximizes the long time scale p-mode effect. In order to study which effect 
introduces larger errors when retrieving the magnetic field vector, inversion 
of simulated profiles are required. 
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Fig. 4 Polarization cross talk induced by the short time solar p-mode effect for 
the three polarization modulations under consideration: scheme A (top panels), 
scheme B (middle panels) and scheme C (bottom panels). Horizontal axis refer to 
the original polarization levels (in percentage of the continuum intensity units). 
Vertical axis indicate the polarization levels after the p-mode effect has been taken 
into account. Left panels are for the linear polarization: y/ Q'^ + . Right panels 
for the circular polarization: V . See text for further details. 



4 Inversion strategy 

4.1 OBSERVED PROFILES 

Prior to the inversion, we must obtain a sufficiently large number of ob- 
served Stokes profiles whose atmospheric parameters are known. To this end 
we have built up a data base of atmospheric parameters with more than 
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20 000 atmospheres. Each of those is represented by a set of parameters 
Xr = [B, 7, i/i, v'^^^, S, S, a]. Each is used to produced a synthetic (that 
simulates observed ones) Stokes vector: I°'''*(Aj) sampled in 300 wavelength 
positions to approximate a continuos function (Equation (3)). This is done 
according to Equation (10). Note that Xj. is used to synthesize the magnetic 
atmosphere: I™s(Aj). To this profile we add a non-magnetic contribution 
P™(Aj) whose shape is known. 

A random time (to) and a modulation matrix {M) must also be chosen. 
According to this and Equation (3) the HMI filter profiles Gi are applied and 
l°^^{Xi,to) is obtained at wavelength position Aj. A linear combination of the 
components of I°'^^(Aj, tg) is produced through Equation (1): Op{Xi,to). The 
next linear combination Op+i{Xi,to + At) at the same wavelength position 

is obtained out of l°^^{Xi,to + At), where these profiles have been shifted 
in wavelength to account for the variation over that time interval At of the 
function ^{t) (see Equation (12) and Figure 2). Note that the only velocity 
that changes in time is the p-mode velocity ^(t). In Equation (12) Vpi^^jj^g^ 
corresponds to v'"'-'^ from Xj- and does not change during the observing 
time. Unless otherwise specified, photon noise is also added to each Op{Xi, to) 
according to Section 3.1. 

In Figure 2, the rms amplitude of the p-mode velocity is about 250 m s~^. 
This was obtained from a two-hour MDI temporal series in the quiet Sun 
using the high resolution field of view. For magnetized regions the amplitude 
of the solar p modes is known to decrease, however. To this end, the rms 
value of the ^{t) function in Figure 2 is scaled with the magnetic field B from 
Xr according to (see Norton, 2000): 



m _ / 250 - 0.155 X B if B< 1000 G\ . _i, , . 
W4*J - 1 128 - 0.033 X B if B> 1000 G'/ ^™ ^ J ^ 

Once all needed Op{Xi), with p = 1, P, are known at wavelength posi- 
tion Aj, the final observed vector I°'''*(Ai) is obtaining by demodulating as in 
Equations (4) - (5). An example of I°'"'(Aj, to) is given in Figure 5 (top four 
panels; dashed line). In circles we also plot the observed profiles, i°^^{Xi,to), 
after the filters Q have been applied. In the lower four panels in Figure 5, 
I°^^{Xi,to) is plotted again. Here we added I°^''(Ai), which arc the observed 
profiles after the effects of the solar p modes have been introduced. Modula- 
tion schemes A (dotted lines), B (dashed lines) and C (dashed-dotted lines) 
have been used. In the three cases, the same to was employed. This plot 
also shows how the solar p modes affect the circular polarization more when 
scheme A is used, and have a stronger effect in the linear polarization profiles 
when scheme B is employed. 



4.2 INVERSION DETAILS 

Once HMI observed profiles I°'^*'(Ai) are available, the next step is to try to 
recover the photospheric parameters Xr that produced them. The underlying 
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Fig. 5 Top 4 panels: Normalized original Stokes profiles I°'"^(Aj) evaluated in 300 

wavelength points (dashed lines). Profiles I°'"'(Ai) after convolution with HMI filter 
profiles (Equation (3); circles). Bottom 4 panels: circles as in the top panels. Dot- 
ted, dashed and daslicd-dotted lines corresponds to the profiles observed by HMI, 

I°'"^(Ai) (see Equation (4)), when modulation schemes A, B, and C, respectively, 
are used. No photon noise has been considered here. 



idea of this inverse process consists in producing a initial synthetic Stokes 
profiles Fy"(Ai) out of an initial set of photosphcric parameters Xq by solving 
the analytical Milne- Eddington solutions to the RTE. The synthetic profiles 
are compared to the observed ones and Xq is modified until the observations 
are matched. The way this works is by building a or merit function as 
follows: 
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2 

wl (21) 

where L is the total number of free parameters {i.e.: seven parameters defin- 
ing Xr), and the indexes k and i run for the four components of the Stokes 
vector and for all wavelengths in which the spectral lino is sampled M =^ 5. 
CTjfc is a measure of the errors in the observations {i.e.: photon noise) and Wk 
is a factor used to ensure that the contribution on the derivatives of with 
K^spcct to the free parameters is not dominated by derivatives with respect 
to Stokes I. To this end, we adopt: 



X 



^ 4 M 

= 7EE 



4°"(A.)-/r(X,A,) 



k=l i=l 



= rrr /^^9n (22) 

'= max{[Ik{X)^]} 

Once is built, we numerically calculate the derivatives of the merit 
function with respect to the free parameters which are then used by a non- 
linear Levenberg-Marquardt least squares fitting algorithm (Press et al. 1986) 
together with a modified Singular Value Decomposition method (Ruiz Cobo 
and Del Toro Iniesta, 1992). The algorithm returns the vector of free pa- 
rameters Xf = [Bf , 7t, tjj[, Vf^^, Sf , Sf, ttf] that minimizes and defines the 
actual {i.e.: real) solar atmosphere: Xr = X/. 

As pointed out, the inversion requires the use of an initial guess model 
Xq. In order to make the simulation as realistic as possible, we choose to fix 
the initial set of parameters (see Table V) that will be kept for all inverted 
profiles. In that way we assume no previous knowledge of the real Xr. Besides 
the initial set of parameters (Xo) wc run a number of extra inversions of 
each Stokes vector hut using randomly generated X's, and only the one that 
retrieves the smallest is kept. 

It is important to mention that our set of free parameters (see Table 
V) does not include the Doppler width, line-to-continuum opacity ratio and 
damping (parameters also used by the Milnc-Eddington model). This strat- 
egy has been adopted to mimic the number of free parameters of a full inver- 
sion method (see for instance Ruiz Cobo and Del Toro Iniesta, 1992) where 
the neglected parameters are determined through the temperature stratifica- 
tion T(r). 



5 Effect of Polarization Modulation 

In this section we will study which polarization scheme allows for a better 
recovery of the physical photospheric parameters Xr and in particular of the 
magnetic field vector. To this end, three simulations (see Section 4 for details) 
were carried out using a different modulation matrix M (see Equations (6) - 
(8)). Each of them includes the inversion of 20 000 individual profiles. Those 
profiles were produced using a set of photospheric parameters (Xr), and a 
initial random time to to take into account the solar p modes (as in Figure 
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Table V Initial set of physical parameters Xo used in the inversion. We only 
assume a previous knowledge in the magnetic field inclination, where 7 < 90° or 
> 90° is chosen depending on the sign in Stokes V lobes. Note that the Doppler 
width, opacity ratio and damping parameters are not considered as free parameters 
(see text for details). 



Physical Parameter 


Symbol 


Value 


Magnetic field strength 


B 


1000 Gauss 


Magnetic field inclination 


7 


60 - 120° 


Magnetic field azimuth 




30° 


Line-of-sight velocity 


yLOS 


2.5 km s"^ 


Source Function 


S 


0.3 


Source Function Gradient 


s 


0.45 


Fining Factor 


a 


0.5 



2; top panel). The p-mode function ^{t) was scaled according to Equation 
(20). At this point, photon noise was not taken into account. The inversion 
of each profile returns a set of inverted atmospheric parameters (Xf) that 
can be compared to the real one (Xr). This allows us to study the error in 
the retrieval of those parameters. 

Figure 6 displays the errors in the retrieval of magnetic field strength 
(top left panel), magnetic field inclination (top right panel), magnetic field 
azimuth (bottom right panel) and linc^-of-sight velocity (bottom left panel). 
The errors are presented as a function of the physical parameter itself. This 
measure is taken by sorting the differences between the real and inverted 
values in increasing order and picking the value that leaves 68 % of the 20 
000 inverted profiles below it. This provides a robust measure of the standard 
deviation in the inferred parameters and will be used hereafter as error. 

In spite of having the longest observing sequence {30 At seconds) , scheme 
C minimizes the errors in the retrieved parameters. This indicates that the 
short term p-mode effect is more important than the long one (see Section 
3.2). In this sense, scheme C ensures that the induced cross talk is minimized 
since orthogonal observations are done with only t = At = four seconds delay 
(see Equations (13) -(14)). 

— Magnetic field strength: scheme A does a better job that scheme B. The 
main reason for this is to be found in the generally smaller amplitudes 
of Stokes Q and U with respect to Stokes V. The same amount of cross 
talk will alter more significantly the linear polarization than the circular, 
thereby affecting in a more severe manner the recovered values for B± 
in scheme B than the retrieved B|| in scheme A. Note that all schemes 
show an increase in the errors when the field strength grows above 1500 
Gauss. This is related to saturation effects: for larger fields, the Zeeman 
splitting is large enough to shift the spectral line out of the scanned 
region, resulting in a loss of information. For small field strengths the 
high increase in the errors, in particular for schemes A and B, is mainly 
due to the larger amplitude of the p modes (Equation (20)). 

— Magnetic field azimuth: as expected, scheme A retrieves smaller errors 
than scheme B. The reason is that the field azimuth angle depends on the 
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Fig. 6 Errors in the retrieval of the different physical parameters as a function of 
the parameter itself. Top left panel: magnetic field strength; Top right panel: mag- 
netic field inclination; Bottom left panel: line-of-sight velocity; Bottom left panel: 
magnetic field azimuth. Different line styles correspond to: modulation scheme A 

solid; Equation (6)), scheme B (dotted; Equation (7)), scheme C (dashed; Equation 

8)). 



ratio ~ Q/U (JefFeries and Mickey, 1991). Since scheme A minimizes the 
p-mode cross talk into the Unear polarization (see Figure 4), it retrieves 
more accurate azimuthal angles. 
— Magnetic field inclination: scheme B gives slightly superior performance 
as compared to scheme A. This has no trivial explanation, since this pa- 
rameter depends on both linear and circular polarization. In addition, all 
three modulation schemes display this kind of behavior in cr-y as a func- 
tion of 7. The error increases towards 7^0°, 180°, while it is minimum 
at 7 ~ 90°. 

Another interesting comparison can be made by plotting the error as a 
function of the filling factor of the magnetic component (a). This is done 
in Figure 7. Errors increase with decreasing a for all physical parameters 
and for all modulation schemes. The reason in this case is trivial. The closer 
a is to 0, the less amount of total observed light comes from the magnetic 
atmosphere I™^ (see Equations (10) -(11)), therefore becoming increasingly 
more diflacult to infer its properties. Hereafter we will focus in the modulation 
schemes A and C, since scheme B can be ruled out due to the larger errors 
that introduces. 
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Fig. 7 Same as in Figure 6 but as a function of the filling factor of the magnetic 
component [a). 



6 Effect of the Photon Noise 

The simulations in Section 5 did not include photon noise. Therefore, the 
next step in our analysis is to consider the simultaneous effect of the solar 
•p modes and the photon noise. To this end we have carried out two new 
simulations where HMI profiles were obtained using modulation schemes A 
and C and taking into account the expected photon noise (sec Section 3.1). 
Errors as a function of the filling factor of the magnetic component (a) are 
presented in Figure 8 in thick lines (solid for scheme A and dashed for scheme 
C). For comparison purposes we plot also the results from schemes A and C 
from Section 5 (where only solar p modes were considered) in thin lines. By 
comparing thin and thick lines in Figure 8 one can see that errors introduced 
by the photon noise are dominant over those induced by the solar p modes. 

An interesting effect appears after introducing the photon noise. As al- 
ready mentioned in Section 5, the errors (when only p modes are considered) 
decrease with decreasing filling factor. This is generally also the case when 
•p modes and photon noise are considered simultaneously. However, for large 
filling factors, a ^ 1, this effect is less pronounced. For some physical pa- 
rameters, like the field strength (Figure 8; top left panel) errors even tend to 
slightly increase. This might be related to the fact that in our simulations, 
those points of the database with large filling factors are intended to repre- 
sent umbral points. The reduced light level in those regions translates into a 
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Fig. 8 Errors in the retrieved physical parameters as a function of the fiUing factor 
of the non-magnetic component. Thin hnes correspond to the results when only p 
modes are considered (same as in Figure 7). Thick lines corresponds to results when 
solar p modes and photon noise (pn) are taken into account. Schemes A and C are 
indicated in solid and dashed lines, respectively. 



smaller signal-to-noise in the observed profiles, thereby increasing the errors 
in the retrieval of the magnetic properties. 



7 Temporal Averages 



In Section 3 we pointed out that at a given wavelength the crosstalk induced 

by the solar p modes changes sign as the oscillation goes from the increasing 
velocity regime to the decreasing part of the wave (sec Figure 3). By averaging 
in time the observed signal, the p-mode crosstalk should be minimized. This 
effect is shown in Figure 9 where we plot in circles the original profiles as if no 
p modes existed: I°^^{Xi). Overploted in dashed lines are Stokes V observed 
during consecutive time intervals: [to, to + 4PM], [to + 4PM, to + 8PM], etc. 
They are computed considering the p modes. These dashed profiles recorded 
during PMAt seconds will be referred to as building block profiles. Since they 
are consecutive in time, the time-dependent behavior of the solar oscillation 
(see Figure 2; top panel) is sampled, each of them contributing differently 
to the induced crosstalk. When all building blocks are averaged, the result- 
ing profile (solid lines in Figure 9) closely resembles that without p-mode 
contribution (circles). 
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Table VI Main properties of the simulations carried out in Section 6. First column 
indicates the modulation scheme used. The asterisk indicates a simulation done 
with the minimum time required to sample the fuU Stokes vector for that particular 
scheme: building block profiles. Second column refers to the amount of total time 
used to average profiles. Note that it is always an integer multiple of the building 
block. Third column is the total number of times that the full Stokes vector I is 
sampled during the time indicated in the second column. Last column shows the 
noise level (ai) in the averaged polarization profiles. This is indicated in units of 
noise level in the building block observation using scheme C. 



Scheme 


Averaging time [secj 


# Stokes I 


(Ti 


A 


880 


10 


0.369 


A 


560 


7 


0.462 


A 


240 


3 


0.707 


A* 


80 


1 


1.22 


(J 


840 


7 


0.377 


C 


600 


5 


0.447 


C 


240 


2 


0.707 


c* 


120 


1 


1 



Figure 9 is obtained from a total averaging time of 480 seconds (roughly 
1.5 oscillatory periods). Note that, as imposed by Equation (6) and Equation 
(8), for modulation scheme A (left panel) the individual observations (dashed 
lines) are much more affected by the solar p modes than for scheme C (right 
panel). This does nothing but to stress the results from Figure 4. Interest- 
ingly, the averaged Stokes V (solid line) arc in both cases very similar. The 
reason for this is that during 480 seconds, scheme A is able to obs(TV(^ Stokes 
V as many as six times: _P = 4 ^ t = 80 seconds. Scheme C is able to record 
Stokes V only four times: P = 6 ^ t = 120 seconds. This allows scheme A 
to provide a better sampling of the solar p-mode oscillations, thus reducing 
more efficiently the induced crosstalk. 

Therefore, it is important to quantify how results compare when using 
different averaging times. In particular, whether scheme A becomes at some 
point more reliable than scheme C. To this end we have carried out several 
simulations using different averaging times with modulation schemes A and 
C. In these simulations photon noise was added to the building block profiles 
according to Section 3.1. The averaged profiles obviously present a larger 
signal-to-noise. Some important properties of the simulations are presented 
in Table VI. 

Figure 10 shows the errors in the physical parameters as a function of the 
time used to average. Note that the temporal evolution of solar structures 
over the averaging time is neglected (except for the solar p modes). In this 
sense the retrieved physical parameters should be considered as an average 
of the physical conditions present in the solar atmosphere during the time 
the data is acquired. Results are presented considering three different solar 
structures: sunspots, plages and network regions. Those definitions arc given 
according to the filling factor of the magnetic atmosphere. Sunspots: a € 
[0.7,1], Plage: a e [0.4,0.6], Network: a e [0.1,0.3]. With the exception of 
the velocity, scheme A (thin lines) becomes more reliable than scheme C 
(thick) using an averaging time as short as 240 seconds. For large averaging 
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Fig. 9 Left panel: example of Stokes V profiles as observed by HMI neglecting the 
effect of the solar p modes (circles) . Dashed lines represent individual observations 
in consecutive time intervals where solar p modes were accounted for using scheme 
A. Each observation requires 5PAt seconds. The average of all of these observations 
is represented by the solid line. Right panel: same but using scheme C. Photon noise 
was not included in these examples. 



Region 


CTB 








[Gauss] 


[deg] 


[deg] [m s"^] 


Sunspot 


4-5 


0.2 


0.2 20 


Plage 


10 


1 


1 30 


Network 


40 - 50 


2 - 3 


5 - 7 60 - 90 



Table VII Approximate errors in the determination of the properties (magnetic 
field vector, velocity) of the magnetic component. This corresponds to a ten minute 
temporally averaged data product. 



Region 


OB 










[Gauss] 


[deg] 


[deg] 


[m s-i] 


Sunspot 


30 


1 


1 


40 


Plage 


80 - 100 


4 


6 - 8 


100 - 200 


Network 


400 - 700 


30 


30 - 50 


> 1000 



Table VIII Approximate errors in the determination of the properties (magnetic 
field vector, velocity) of the magnetic component. This corresponds to a 120 second 
data product (building block profiles) using Mod C. 



times, both observing schemes become almost identical. The HMI scientific 
goal is to provide the full magnetic field vector with a cadence of ten minutes. 

At this level the errors for sunspot regions could be as little as (Tb ^ 4 Gauss, 
cr^ ~ 0.2°, ~ 0.2° and (7„ ~ 20 m s'^ (see Tables VII and VIII for details). 



8 Conclusions 

The HMI instrument intends to provide accurate and continuous full-disk 
measurements of the solar magnetic field vector in the photosphere. We have 
identified and carefully studied the main sources of errors in HMI: limited 
spectral resolution, solar p modes (that affect our data as a consequence of the 
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Fig. 10 Errors in the determination of the different physical parameters as a func- 
tion of the averaging time. Thin lines indicate results for modulation scheme A 
whereas thick ones are for scheme C. Different line styles correspond to differ- 
ent observed solar regions; sunspots (solid), plage (dotted) and network regions 
(dashed). 



non-simultaneity in the polarization measurement^ and photon noise. The 
most critical source of error turns out to be the photon noise. Increasing the 
signal-to-noise ratio of the observations is achieved with longer integration 
times. This is not trivial. Increasing the exposure time per measurement 
may lead to saturation effects as we reach the limits of the CCD's linear 
behavior. The only possibility left is to combine several measurements. We 
find that averaging data over ten minutes increases the signal-to-noise by a 
factor of two to three (see Table VI) plus effectively removes the effect of 
the solar p modes. This leads to very precise determinations of the time- 
averaged magnetic field vector, with an accuracy better than five Gauss in 
field strength and 0.2 degrees in inclination and azimuth. 

A faster and better correction for the p modes can be achieved in different 
ways. On the one hand, the observations can be averaged over time using 
more sophisticated methods. In Section 7 a simple (unweighted) average was 
used. More sophisticated ways {e.g.: using appropriate weights) would lead 
to a faster smearing of the solar p modes. On the other hand, it is possible to 
interpolate the data so that we combine polarization measurements taken at 



^ Note that non-simultaneous observations will not be affected by seeing induced 
cross talk since HMI is a space-borne instrument 
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the same time (see Equation (14)). This might correct the effects of the solar 
p modes even without time averaging. A last possibility would be to apply a 
subsonic filter to the raw {i.e.: monochormatic and not demodulated) data 
(Rutten, Wijn, and Siitterlin, 2004). Although these possibilities are worth 
investigating, they will not significantly change our results since the main 
source of errors is the photon noise. 

A critical point we have not mentioned concerns the contribution of the 
non-magnetic atmosphere I™ (see Equation (10)). In all our simulations 
we have considered I"™ as known (with the exception of the filling factor 
a). This is a best case scenario. In reality, the I""^ contribution can arise 
because of scattered light from the optical elements of the instrument, or 
due to the presence of a real non-magnetic atmospheric component within 
the resolution element of the observations. 

In the first case (instrumental scattered light), it might be possible to 
characterize I"™ with laboratory tests. In the second case, its origin would 
be solar, and therefore we should have considered new free parameters in our 
inversion to describe I"™, such as: vji;^^, Snm, Snm- We have neglected this 
in order to keep the number of free parameters (listed in Table IV) within 
reasonable limits and below the number of data points (20). A necessary ex- 
tension of this work should address this point carefully by both (a) including 
the new free parameters and, (h ) by using a different I'™ in the synthesis and 
the inversion so that we can quantify the errors introduced by a poorly known 
Obviously the larger a, the smaller these errors will be. Therefore, for 
sunspot regions (Figure 10 and Tables VII and VIII) the results presented in 
this work should be robust. 

Another very important issue concerns the speed of the inversion process 
(see Section 4) . The employed method is a traditional non- linear least-squares 
fitting algorithm. This method will not ultimately be used for the inversion 
of HMI data. Our simulations with the analysis of 20 000 profiles took an 
average of 24 hours to complete. HMI requires the processing of millions of 
profiles per minute. Even those methods, considered to be fast, such as Prin- 
cipal Component Analysis (Rees et ai, 2000) or Artificial Neural Networks 
(Socas-Navarro, 2003), will be challenged by HMI requirements. This will be 
the subject of a future investigation. 
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